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We report results on time-dependent CP asymmetries in B — > D'^tt^ decays based on a data 
sample containing 657 x 10 6 BB pairs collected with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider at the T(4S) resonance. We use a partial reconstruction technique, wherein 
signal B — > D*^^ events are identified using information only from the fast pion from the B decay 
and the slow pion from the subsequent decay of the D* T , where the former (latter) corresponds to 
D* + (D*~) final states. We obtain CP violation parameters S + = +0.061 ±0.018(stat) ±0.012(syst) 
and S~ = +0.031 ± 0.019(stat) ± 0.015(syst). 

PACS numbers: 11.30.Er: 14.40.Nd 



In the Standard Model (SM), CP violation occurs due to the presence of a complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [j. Precision measurements of the parameters of the CKM matrix are important to inves- 
tigate new sources of CP violation. The study of the time-dependent decay rates of B°(B a ) — > D*^^ provides 
a theoretically clean method for extracting sin(20i + 03) 0, where 0i and 03 arc angles of the CKM Unitarity 
Triangle as defined in 0. As shown in Fig. [TJ these decays can be mediated by both Cabibbo-favored (CF) and 
doubly-Cabibbo-supprcsscd (DCS) diagrams, whose amplitudes are proportional to V* b V u d and V* b V c d, respectively, 
where Vij are the CKM matrix elements and have a relative weak phase difference 03 . 

The time-dependent decay rates are given by [I| 



P(B° -> D*^) 
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|At|/r B o 
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[l + Ccos(AmAt) - S' ± sin(AmAi)] , 

n*±^=F^ = 1 e -\M\/r B o 



P(B° -> D^ttF) 



x [l ± C cos( AmAt) + S ± sin( Am At)] (1) 

Here At is the difference between the time of the decay and the time that the flavor of the B meson is tagged by 
the associated B meson, Tgo is the average neutral B meson lifetime, Am is the B°-B° mixing parameter, and 
C = (l — R 2 ) j (l + i? 2 ), where R is the ratio of the magnitudes of the DCS and CF amplitudes (we assume their 
magnitudes to be the same for B° and B° decays). The CP violation parameters for D*ir are given by 

± _ -2flsin(2^+fo±g) ^ 

b ~ (TTW) ' (2j 

where 6 is the strong phase difference between the CF and DCS amplitudes. 

d / -c 

b * < c b £ . u 

B° D " B° 7t + 

d > d d > d 

FIG. 1: Diagrams for B° -»■ D*-% + (left) and B° -> D* _ 7r + (right). Those for B° -> D'+tt' and B° D* + -k~ can be 
obtained by charge conjugation. 

Since the predicted value of R is small, ~0.02 Q, we neglect terms of O (R 2 ) (and hence take C = 1). The amount 
of CP violation in D*tt decays, which is proportional to R, is expected to be small, and hence, a large data sample is 
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needed in order to obtain sufficient sensitivity. To increase statistics, we empioy a partiai reconstruction technique [f|, 
wherein signal is distinguished from background on the basis of kinematics of the 'fast' pion (itf) from the decay 
B —> D*irf, and the 'slow' pion (w s ) from the subsequent decay of D* —> Dtt s ; thus the D meson is not reconstructed 
at all. 

Previous analyses have been reported by Belle @, as well as by BaBar [{|. This study uses a data sample of 
605 fb -1 containing 657 x 10 6 BB events. The data sample is about twice the size of the dataset used in the previous 
Belle analysis Q and supersedes the previous study. 

The data were collected with the Belle detector [lj| at the KEKB collider operating near the T(4S) resonance. 
The Belle detector is a large-solid- angle magnetic spectrometer that consists of a silicon vertex detector (SVD), a 50- 
layer central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel- like arrangement 
of time-of-fhght scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoidal coil that provides a 1.5 T magnetic field. An iron flux- return located 
outside of the coil is instrumented to detect K L mesons and to identify muons (KLM). A sample containing 152 x 10 6 
BB pairs was collected with a 2.0 cm radius beampipe and a 3-layer silicon vertex detector (SVD1), while a sample 
of 505 x 10 6 BB pairs was collected with a 1.5 cm radius beampipe, a 4-layer silicon vertex detector (SVD2), and a 
small-cell inner drift chamber (l2j . 

The "signal side" B, decaying to D* + irJ , D* + — > D a irf (or charge conjugate), is reconstructed using pairs of 
oppositely charged pions. Since the pion originating from the B has a higher momentum in the T(4S*) cm. frame 
than that originating from the D* , the former (latter) is referred to as the fast (slow) pion. All momenta and energies 
in this paper are calculated in the T(4S') center-of-mass (cm.) frame, unless otherwise stated. Fast pion candidates 
are required to have a radial (longitudinal) impact parameter dr < 0.1 cm (\dz\ < 2.0 cm) and to have associated hits 
in the SVD. We reject leptons and kaons based on information from the CDC, TOF and ACC. A requirement is made 
on the fast pion momentum, 1.93GcV/c < pf < 2.50GcV/c. Soft pion candidates are required to have momenta 
in the range 0.05GeV/c < p s < 0.30GcV/c. No particle identification requirement is applied for these pions. We 
impose only a loose requirement that they originate from the run-dependent interaction point (IP) profile. The IP 
has a z ~ 4mm along the beam direction (z), and a x ~ 100 /im and <j y ~ 10 /im in the plane perpendicular to the 
beam direction. 

For any given tt f from a signal B decay, the energy of the D* may be known through energy conservation, Ed* = 
Eb — E 7Tf , where Eb = s/s/2 at the T(4S). The magnitude of the momentum is then \pd* \ = \J E' Dr — rri D , . Because 

the B meson is slow in the cm. frame, its momentum \ps\ = J Eg — m 2 BQ « 0.3 GeV/c is small relative to the Hf 
and D* momenta. It follows from momentum conservation 

PD' = Pb - Pit i (3) 

that the direction of the D* momentum can be approximated as the direction opposite to p n . . This approximate D* 
four-momentum is denoted as the "partially reconstructed" D* . We define a quantity ps = \pn f \ ~ which for 

signal decays satisfies \ps\ < \pb\, as can be seen by examining Eq. ([3]). 

We then examine the soft pion after boosting it into the partially reconstructed D* frame; in the true D* rest 
frame, the soft pion is monoenergetic and its momentum has an angular distribution characteristic of a pseudoscalar 
to pseudoscalar-vcctor transition, oc cos 2 9 where 6 is taken relative to the boost axis. In the partially reconstructed 
frame, the momentum will have a limited spread. We study the components parallel and perpendicular to the boost 
axis, denoted p|| and p±, respectively. 

We use the three kinematic variables ps, p» and p± to distinguish between signal and background. Background 
events are separated into three categories: D* =F p ± , which is kinematically similar to the signal; correlated background, 
in which the soft pion originates from the decay of a D* that in turn originates from the decay of the same B as the 
fast pion candidate, excluding Z?* T 7r ± and D* =F p ± decays {e.g., B -> D**n, B -» D*a\, B — > D*lv); and uncorrelated 
background, which includes all other background sources (e.g., continuum processes, B — > Dir). The distributions 
of the kinematic variables for signal and background categories are determined from a large sample of Monte-Carlo 
(MC) generated data corresponding to three times the integrated luminosity of our data sample. 

We retain candidates that satisfy -0.10 GcV/c < py < 0.07 GeV/c, -0.60 GcV/c < ps < 0.50 GcV/c and p± < 0.05 
GeV/c. In the cases where more than one candidate satisfies these criteria, we select the one with the largest value 
of S 7rf7! - g , where <5 7r/7rs is the angle between the fast pion direction and the soft pion direction in the T(4S') cm. 
frame. The signal region is defined as: -0.40GeV/c < ps < 0.40GeV/c, -0.05GeV/c < p\\ < -O.OlGeV/c or 
0.01 GeV/c < p|| < 0.04GeV/c and p ± < 0.05 GeV/c. 

The determination of the flavor of the B meson opposite to the signal side B, which we refer as the tag-side B, 
is essential for this measurement. In order to tag the flavor of the associated B meson, we require the presence of 
a high-momentum lepton (I) in the event. This helps reduce background from continuum e + e~ — > qq (q — u, d, s, c) 
processes. Tagging lepton candidates are required to be positively identified cither as electrons, on the basis of 
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information from the CDC, ECL and ACC, or as muons, on the basis of information from the CDC and the KLM. 
They are required to have momenta in the range 1.1 GeV/c < pi < 2.3 GeV/c, and to have an angle with the fast 
pion candidate that satisfies cos(5 7r/ ; > —0.75 in the T(4S) cm. frame. These requirements reduce to a negligible 
level (0.7%) the contribution of leptons produced from semileptonic decays of the unreconstructed D mesons in the 
B -> _D* =F 7r ± decay chain. 

Vertexing requirements identical to those for the fast pion are applied to the lepton candidate in order to obtain 
an accurate vertex position. To further suppress the remaining small continuum background, we impose a loose 
requirement on the ratio of the second to zeroth Fox- Wolfram [13J moments, i?2 < 0.6. 

Event-by-event signal and background fractions arc determined from binned maximum likelihood fits to the two- 
dimensional distributions of p$ and pu . The results of these fits, projected onto each of the two variables, are shown 
in Fig. [5J and summarized in Table Q] We obtain a purity of 59.0 ± 0.4% in the signal region, where purity is defined 
as the ratio of the signal to total yields. 




FIG. 2: Results of the fits to D*n candidates projected onto pu (left) and ps (right) in the signal region of the two kinematic 
variables. The contributions are: D*ir (open), D* p (green slanted lines), correlated background (magenta crossed lines) and 
uncorrelated background(shaded blue). Data are shown as points with error bars. 



TABLE I: Summary of the yields in the signal region 



D*tv 50196 ±286 

D*p 10232 ± 150 

Correlated background 10425 ± 135 
Uncorrelated background 14193 ± 128 



At the KEKB asymmetric-energy e + e~ (3.5 GeV on 8 GeV) collider, operating at the T(4S) resonance (E c . m = 
10.58 GeV), the T(4S) is produced with a Lorentz boost of /?7 = 0.425, almost along the electron beamline (z). In 
the T(45) cm, B° and B° mesons are approximately at rest. Hence the proper timc-differcncc (At) between the 
signal side vertex (z s ; g ) and the tag-side vertex (ztag) is obtained from the fast pion on the signal side and the tagging 
lepton. The variable At is defined as: 

At « (z sig - z tag )/^7c- (4) 

z s ig is obtained from the intersection of the fast pion's track and the IP, and Ztag is obtained from the intersection of 
the tagging lepton's track and the IP. 

To measure the CP violation parameters, we perform a simultaneous unbinned fit to four samples: two are of 
same-flavor (SF) events, namely Tt + l + , , in which the fast pion and the tagging lepton have the same charge, and 
the other two are of opposite-flavor (OF) events, namely it + 1~ , 7r~Z + , in which the fast pion and the tagging lepton 
have opposite charge. We minimize the quantity — ln£ = — ln£i, where 

£i fD*7rPD*7r fD*pPD*p ~T" ,/unco-funco ~t~ /corr-fcorr- (5) 

Here, f x stands for the event-by-event fraction from source x and is obtained from the fits to the kinematic variables, 
and P denotes the probability density functions (PDFs) for signal and backgrounds, which contain an underlying 
physics PDF with experimental effects taken into account. The convolution of the physics PDF with experimental 
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effects will be described later. For D*tt and D*p, the PDF is given by Eq. (Q}, whereas for D*p the terms are 
effective parameters averaged over the helicity states [3] and are constrained to be zero. The PDF for correlated 
background contains a term for neutral B decays (given by Eq. |T|) with S ,± = 0), and a term for charged B 
decays (for which the PDF is 2r 1 - e~l At l/ r -B+ , where t b + is the lifetime of the charged B meson). The PDF for 

uncorrelated background also contains neutral and charged B components, with the remainder from continuum 
e + e~ — > qq (q = u,d,s,c) processes. The continuum PDF is modeled with two components: one with negligible 
lifetime, and the other with a finite lifetime, which takes into account the dependence of average lifetime of the charm 
contribution in the continuum (close to the average D meson lifetime). 

The parameters in P U nco and P C orr are obtained from separate simultaneous fits to OF and SF candidates in 
the respective sideband regions, defined later. Since there is no CP violation in background, the corresponding 
parameters are fixed to zero in these fits. The fit is further simplified by fixing the biases in Az to zero (discussed 
later in detail). MC simulation studies demonstrate that varying or fixing these biases to zero does not affect the 
background parameters. 

To measure the uncorrelated background shape, we use events in a sideband region, — O.lOGeV/c < pn < 
-0.07 GeV/c or 0.01 GcV/c <py < 0.04GcV/c, -0.60 GeV/c < p s < 0.50 GeV/c and 0.08 GeV/c < p ± < O.lOGcV/c, 
which is populated mostly by uncorrelated background (~ 90%). To determine the correlated background parame- 
ters, we use events in a sideband region, — O.lOGeV/c < p» < — 0.07GeV/c, — 0.60GeV/c < ps < O.OOGeV/c and 
0.00 GcV/c < p± < 0.05 GeV/c. This sideband region is dominated by both correlated and uncorrelated backgrounds 
and has a very small amount of D*ir signal and D*p background. The uncorrelated background parameters are 
fixed to the values obtained in the previous fit. Figure [3] shows p± distributions for signal and various background 
components in MC simulations, corresponding to about three times the size of the data. 




FIG. 3: p± distributions for various MC simulations, showing the following contributions: D*n (red solid line), D* p (dotted 
green line), correlated background (dashed magenta line), and uncorrelated background (blue solid line). 

The PDF for signal and background in Eq. ([5]) must be convolved with the corresponding Az resolution functions 
related to the kinematic smearing (TZk), detector resolution (T^dct), an d asymmetry in Az from non-primary tracks 
(72-np)- The resolution function related to kinematic smearing is due to the fact that we use the approximation of 
Eq. ((¥]). The detector resolution and smearing due to the asymmetry in Az from non-primary tracks are described in 
detail elsewhere [1]. 

To account for mistagging, the PDFs in Eq. ([5]) arc divided into two components 

P(l T ,irf) = (1 - w T )P(B°/B° -> D'^tt*) 

+w±P(B°/B° -> ZT^), (6) 

where w + and w~ are the wrong-tag fractions, defined as the probabilities to incorrectly measure the flavor of tagged 
B° and B° mesons, respectively, and are determined from the data as free parameters in the fit for S . 

The time difference At is related to the measured quantity Az as described in Eq. (j4|) , with an additional term due 
to possible offsets in the mean value of Az, 

At^At + e At ~ (Az + e A z) /he. (7) 

It is essential to allow non-zero values of ca* since a small bias can mimic the effect of CP violation: 

cos(AmA^) -> cos(AmAi) - Ame A t sin(AmAt) (8) 

A bias as small as e Az ~ 1 ^m can lead to sine-like terms as large as 0.01, comparable to the expected size of the CP 
violation effect. Because both vertex positions are obtained from single tracks, the partial reconstruction analysis is 
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more susceptible than other Belle CP violation analyses to such biases. We allow separate offsets for Az for each 
combination of and I charges. Thus we have eight offsets in total, four for each data sample, SVD1 and SVD2. 

To extract the CP violationparameters we fix tb° and Am at their world average values (rgo = 1.530 ± 0.009 ps 
and Am = 0.507 ±0.005 ps -1 Q), and fit with S + , S~ , two wrong tag fractions, and eight offsets as free parameters. 
We obtain S + = +0.061 ± 0.018 and S~ = +0.031 ± 0.019, where the errors are statistical only. The wrong tag 
fractions are = (5.3 ± 0.3)% and w+ = (5.2 ± 0.3)%. All floating offsets arc consistent with zero except for one 
of the OF combinations (717 = ir~ , I = £ + ) in the SVD1 sample. The results are shown in Fig. |4] Using large MC 
samples generated with non-zero and zero S values, we do not find any significant bias in the procedure. 

To further illustrate the CP violation effect, we define asymmetries in the same flavor events {A ) and in the 
opposite flavor events (.A OF ), as 

N v - l -(Az)-N n+l+ (Az) 
' N 7r - l -(Az)+N„ +l+ (A z )' 
a of = N v+l -(Az)-N„- l+ {Az) 

' N n+l -(Az)+N w - l+ (A z y w 

where the TV values denote the number of events for each combination of / and I charge. These are shown in Fig. [5] 




FIG. 4: Az distributions for four flavor- charge combinations: ty~1~ (top left) , ty~1 + (top right), ir + l~ (bottom left), and 
ir + l + (bottom right). The fit result (solid blue line) is superimposed on the data (solid points with error-bars). The signal and 
background components are shown as the solid red and dotted black curves, respectively. 




FIG. 5: Results of the fit to obtain S + and S , shown as asymmetries in the SF events (left) and OF events (right). The fit 
results (solid blue lines) are superimposed on the data. 

This analysis is very sensitive to the vertexing bias. Hence, we include Az offsets in the fits to account for this bias. 
In order to estimate the error due to these offsets, we perform fits to obtain values with and without offsets using 
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an ensemble of 100 generated D*n signal samples, and use the difference between the two results as the systematic 
error. We obtain negligible contribution to the systematic errors when we float Az offsets in the background PDF. 

Other sources of systematic error are the resolution functions, 7?.^, 7Z-dct and !Z np , uncorrelated and correlated 
backgrounds and physics parameters, Am, t b o and t b + that are fixed in the fit to extract 5 ,± . The parameters of 
the resolution functions and backgrounds are varied by ±1(7 (with Am and t b o fixed), respectively, where cr's are 
the corresponding errors of the parameters and the difference is assigned as systematic error. We vary the physics 
parameters by ±1(7, where a is the error of the corresponding PDG values, and we then use the difference between the 
values thus obtained and the default values as the systematic error. When the fit is performed floating values, 
along with t b o and Am, we obtain: S+ = +0.055 ± 0.018 and S~ = +0.039 ± 0.019, t B o = 1.550 ± 0.008 ps and 
Am = 0.473 + 0.004 ps -1 , where the errors are statistical only. The deviations from the nominal fit (0.06, 0.08) are 
close to the systematic errors assigned for the physics parameters (Table Hij) . The difference between the values 
obtained floating both Am and t b o parameters and the default value is also added to the systematic error estimation. 
In the fits to extract 5' ± , S ± D' P and S ± COII are set to zero. For the systematic error due to these parameters, the fit 
is performed with these values set to ±0.05 and the difference between the 5 ± value thus obtained and the default 
value is assigned as the systematic error. 

We use a triple Gaussian to model the detector resolution (i?dct) function. We consider the systematic uncertainty 
due to the lack of knowledge of the exact functional form of the resolution model. When the resolution models are 
varied, we obtain shifts as large as 0.006 for S + . This is conservatively assigned as the systematic error due to this 
source. 

We obtain a vertexing systematic error of 0.003 for 5 ± . Additional systematic errors result from varying the number 
of bins for the kinematic variables, p$ and p\\ in the yield fit. 

The systematic errors are summarized in Table [TXJ The total systematic error is obtained by adding the above 
terms in quadrature. 



TABLE II: Summary of possible sources of systematic error 



Systematic error source 


S+ 


s~ 


A 2 offset 


0.002 


0.003 


IZk parameters 


0.002 


0.003 


TZdet parameters 


0.002 


0.002 


1Z np parameters 


0.004 


0.004 


Background parameters 


0.001 


0.001 


Physics parameters 


0.006 


0.009 


Floating t b o and Am 


0.006 


0.008 


Yield fit 


0.003 


0.005 


Resolution model 


0.006 


0.002 


IP constraint 


0.003 


0.003 


Total systematic error 


0.012 


0.015 



In conclusion, we have measured CP violation parameters that depend on <p^ using the time-dependent decay rates 
of B a — > D* =F 7r ± with a data sample containing 657 x 10 6 BB events. We determine the CP violation parameters 
S ± to be 

S + = +0.061 + 0.018 + 0.012, 

S~ = +0.031 + 0.019 + 0.015, (10) 

where the first errors are statistical and the second errors are systematic. We can also express the results as parameters 
a, c, defined as: 

a = + S~)/2, 

c = -(S + -S~)/2. (11) 

Our results thus become: 

a = -0.046+0.011 + 0.015, 

c = -0.015 + 0.011 + 0.015. (12) 

The deviation of a from zero is a measure of the amount of CP violation. We obtain a significance of 2.5a on the 
CP violation parameter, a. Our measurement is consistent with the world average value and significantly improves 
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the precision of previous measurements reported by Belle [7|, |8[ as well as by BaBar [9( and supersedes our earlier 
result [|. 
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